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Abstract 

We consider the problem of finding a lowest cost dominating set in a given disk graph containing 
n disks. The problem has been extensively studied on subclasses of disk graphs, yet the best known 
approximation for disk graphs has remained O(logn) - a bound that is asymptotically no better than 
the general case. We improve the status quo in two ways: for the unweighted case, we show how to 
obtain a PTAS using the framework recently proposed (independently) by Mustafa and Ray [16] and by 
Chan and Har-Peled [4]; for the weighted case where each input disk has an associated rational weight 
with the objective of finding a minimum cost dominating set, we give a randomized algorithm that 
obtains a dominating set whose weig ht is within a factor 2^^^'°^* of a minimum cost solution, with high 
probability - the technique follows the framework proposed recently by Varadarajan [19) . 

1 Introduction 

For a set 2? of rt disks in the Euclidean plane, define an intersection graph, G = {V,E), thus: = P; 
{u, v} E E ■iF^ disk (u) n disk {v) ^0. G is called a disk graph; it is a unit disk graph when the disk radii are 
identical. 

Given a graph the minimum dominating set (MDS) problem is to find a smallest subset V C_V such that 
every vertex is either in V or is adjacent to a vertex in V . On general graphs, the problem is (1 — e) Inn 
hard to approximate for any e > under standard complexity theoretic assumptions |10[ [5] , while a greedy 
algorithm yields an O(logn) approximation [20j . 

Nevertheless, better approximations are possible for restricted domains. For example, the problem admits 
a polynomial-time approximation scheme (PTAS) for unit disk graphs and growth-bounded graphs [131 117j . 
The problem is NP-hard on these domains [B]. However, for the disk graph case, o(logn) approximations 
have remained elusive - perhaps, in part, because known techniques for unit disk graphs and solutions to 
other problems on disk graphs have either relied on packing properties |13[|17i l8l[3]. or when packing property 
does not hold, as in the minimum weighted dominating set on unit disk graphs, the fact that disk radii are 
uniform [T1[TS]. Erlebach and van Lceuwen recently studied the dominating set problem on fat objects, e.g., 
disk graphs, [5]. They note that existing techniques for disk graphs do not seem sufficient to solve MDS [S]; 
they also give an 0(l)-approximation for fat objects of bounded ply. 

In their recent break-through papers, Chan and Har-Peled [4], and Mustafa and Ray pj] independently 
showed how a simple local search algorithm on certain geometric graphs yields a PTAS for some problems; 
Chan and Har-Peled [1] show local search yields a PTAS for maximum independent set problem on admissible 
objects, while Mustafa and Ray [TB] show local search yields a PTAS for the minimum hitting set problem 
given a collection of points and half-spaces in K^, and also for points and admissible regions in R^. They both 
use the planar separator theorem to relate the cost of the local search solution with the optimum solution. 
In the framework, at the crux lies the analysis of a certain graph whose vertices are objects found by local 
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search and ones that belong to an optimum solution, and whose edges (which are only between the two kinds 
of vertices) satisfy a property relating the two solutions. They show that there exists such a graph which is 
also planar. Mustafa and Ray |16j refer to the existence of such a planar graph as the locality condition. 

Results: Our first result is a PTAS for the minimum dominating set problem for disk graphs via a local 
search algorithm, as in pi 116). Our analysis also uses the framework introduced by these two papers. Our 
main new contribution is to show the existence of a planar graph satisfying the locality condition. This 
graph turns out to be the dual of a weighted Voronoi diagram in the plane. 

The minimum dominating set problem for disk graphs can be reduced to the problem of hitting half- 
spaces in with the smallest number of a given set of points. That is, given the set T) of disks that form 
the input to the MDS problem, we can easily compute a map tt from I? to a set of points in M^, and a map 
h from 2? to a set of half-spaces in M^, with the following property: Two disks di and d2 from T) intersect 
if and only if TT{di) lies in h{d2). Thus we can efficiently reduce the MDS problem for disks to a hitting set 
problem for points and half-spaces in M.^. While there is a PTAS for the hitting set problem in R"^, as shown 
by [11], there is none known for R^. It is not hard to see that a local search such as the one in [TB] docs not 
yield a PTAS in M^. 

Rather than reduce to a hitting set problem, we are able to establish the locality condition by staying in 
the plane itself. In fact, the graph for the locality condition is the dual of the weighted Voronoi diagram of 
the centers of the disks in the local search solution and the optimal solution, where the weights arc the radii 
of the disk. This can be seen as generalizing the situation considered by [TB] for the hitting set problem with 
points and disks in the plane. In that case, the graph for the locality condition is the Delaunay triangulation, 
which is the dual of the unweighted Voronoi diagram. 

For the case when the disks are weighted, we give the first o(log n) approximation algorithm; we give a 
20(iog n) approximation algorithn0. This result is based on the framework recently introduced by Varadara- 
jan for the weighted geometric set cover problem |19| . Our contribution here is to observe that the framework 
is applicable to our dominating set problem as well; the weighted Voronoi diagram is the key to this result 
also. 

We assume that the inputs for both problems satisfy non-degeneracy assumptions - no three disk centers 
on a line and no four disks tangent to a circle. This is without loss of generality, as these conditions can 
be enforced by simple perturbations. In Section [2l we present our PTAS for the unweighted dominating set 
problem, and in Section [3] our algorithm for weighted dominating set. 

2 The Unweighted Case: PTAS via Local Search 

In this section, we give our PTAS for minimum dominating set for disk graphs. Here, we are given a disk 
graph with a set T) oi n disks in the Euclidean plane, and we arc interested in computing a minimum 
cardinality dominating set of the disk graph. The algorithm is given in Section 12.11 and the analysis of the 
approximation ratio is given in Section 12.21 

2.1 The Algorithm 

Local Search. Call a subset of disks, BCD, fe-locally optimal if one cannot obtain a smaller dominating 
set by removing a subset X C B of size at most b from B and replacing that with a subset of size at most 
|X| — 1 from 'D\B. Our algorithm will compute a &- locally optimal set of disks for 6 = ^ where c > is a 
large enough constant. Our algorithm begins with an arbitrary feasible set of disks and proceeds by making 
small local exchanges of size b = O(^), for a given e > 0. We stop when no further local improvements are 
possible. 

Suppose that the solution returned is B. Finally, for reasons apparent in the analysis, we check to see if 
for any disk u G B there is a disk v G V such that u is completely contained in v G 'D\ B. If such a disk 

^log* n is the fewest number of iterated "logarithms" appUed to n to yield a constant. 
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exists, then simply replace u with the largest such disk v. We return this as our final solution and call it B. 
Our replacement step ensures that there is no disk in B that is properly contained in some other disk in T). 

Running Time. We will now show that the running time can be bounded by a polynomial in n. The 
number of swaps that the local search algorithm will make is at most n, because there are n disks and each 
swap strictly decreases the number of disks in the solution. For each swap, we need to check every subset of 
disks of size at most h which can be done in time 0(n''). Recall that b is only a function of e, and thus we 
can have the exponential dependence on h in the running time. 

So we will make at most n swaps, each of which takes time 0(n''). Clearly, the last step (where we check 
to see if a disk is contained within another) can be done in time polynomial in n, and therefore the entire 
running time of the algorithm is efficient with respect to n. 

2.2 Approximation Ratio 

We will show that our algorithm is a PTAS, thus proving the following theorem: 

Theorem 1. For any e > 0, there exists a polynomial time algorithm for the minimum dominating set 
problem on disk graphs that returns a solution whose cost is at most (1 + e)OPT where OPT is the cost of 
an optimal solution. 

Let R be the disks in an optimal solution; we may assume no disk in R is properly contained in any 
other disk in T>. Thus, no disk in i? U i? is properly contained in any other disk oi RU B. Note that by the 
definition of PTAS, we need to show that \B\ < (1 + e) • \R\. We will refer to R as the set of red disks and B 
as the set of blue disks. Without loss of generality, we will assume that i? H B = 0, i.e. there is no disk that 
is both red and blue. For a disk u G 2?, we say a disk v € RU B is a dominator of u \i u and v intersect. 
Similarly, we also say that v dominates u. 

Wc must show the existence of an appropriate planar graph which relates the disks in R with the disks 
in B. Here, we state the locality condition as per Mustafa and Ray |16] : 

Lemma 2 (Locality Condition). There exists a planar graph with vertex set RUB, such that for every 
d £ T>, there is a disk u from amongst the red dominators of d and a disk v amongst the blue dominators of 
d such that {u, v} is an edge in the graph. 

Section is devoted to a proof of Lemma We then describe the argument (from JB]) that uses the 
lemma to show that \B\ < (1 + e)|i?|. 

2.3 Establishing the Locahty Condition 

This section is devoted to the proof of Lemma [2l that is, the construction of an appropriate planar graph 
which satisfies the locality condition. 

Weighted Voronoi Diagram. We will be using a generalization of Voronoi diagrams called a weighted 
Voronoi Diagram (WVD). Instead of defining cells with respect to a set of points, we will be defining cells 
with respect to red and blue disks. In order to do this generalization for disks, we must define the distance 
between a point in the plane and a disk. 

Let u be a disk and let a; be a point in the plane. Wc define wvd(a:;, u) — d(x, c^) — r^ where is the 
center of u, r„ is the radius of u, and d{x, c„) is the Euclidean distance between x and c„. Intuitively, for a 
point X, 'W'vd(x,u) is the Euclidean distance from x to the boundary of u; the distance to a disk is negative 
for points that arc strictly inside the disk. Alternatively, if a; ^ u, then wvd(a;, u) is the amount wc would 
need to increase the radius of u so that x lies on the boundary of u; if x € u, then wvd(x, u) is the negative 
of the amount we would need to decrease the radius of u so that x lies on the boundary of u. See Figure [1] 
for an illustration. 
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For a disk u in any collection of disks, let cell(it) be the set of points x in the plane such that wvd(x, u) < 
wvd(x, v),u ^ V. The cells of all the disks in the collection induce a decomposition of the plane, and this is 
the WVD. This is just the standard weighted Voronoi diagram of the centers of the disks, where the weight 
of the center of a disk is simply the radius of the disk P] . 

Consider the WVD of the disks in RUB. First, we will show that for every u G RUB, u has a non-empty 
cell in the WVD. That is, there is some point in the plane that is closer to u than it is to any other red or 
blue disk. 

Lemma 3. In the weighted Voronoi diagram of the union of red and blue disks, the cell of every disk u is 
nonempty. Moreover, Cu (the center of u) belongs only to cell{u). 

Proof. We will show that c„ is only in cell(u). Suppose for the sake of contradiction that c„ € cell(w) such 
that u ^ V. This means that wvd(cii,u) < wvd(cu,it) = d{cu,Cu) — ^ -^r^. So, — r„ > wvd(cti,v) = 
d{cu,Cu) — ry r^ > d{cu,Cu) + ru- This implies that u is contained in v, and since the two disks are not 
the same, the containment is proper. But this is a contradiction, since no disk in RU B contains another 
such disk. □ 



The Graph. Any cell in the WVD of i? U i? is star-shaped with respect to the center of the corresponding 
disk. That is, for every point y e cell(u), the segment c„?/ is contained within cell(M). 

The graph for the locality condition is simply the dual of the WVD oi RU B. That is, for each cell in 
the WVD there is a vertex, and there is an edge between two vertices if and only if their corresponding cells 
share a boundary in the diagram (that is, if and only if there is a point in the plane equidistant from the 
two disks). The graph is planar ~ exploiting the fact that the cells are star-shaped, the edges can easily be 
drawn so that no two edges intersect [2]. 

Corollary 4. The dual of the power diagram of RU B is a planar graph. 

Because every red and blue disk has a nonempty cell in the WVD, every such disk will also have a 
corresponding vertex in our planar graph. We are now ready to show that for each d S there is a disk u 
from amongst the red dominators of d and a disk v amongst the blue dominators of d such that cell(M) and 
cell(i>) share a boundary in the WVD. This would then imply that their corresponding vertices in the graph 
share an edge, completing the proof of Lemma [2l For simplicity, if there is an edge connecting the vertex 
corresponding to cell(M) and the vertex corresponding to cell(v), then we will simply say there is an edge 
connecting u and v. 

Lemma 5. In the dual graph of the weighted Voronoi diagram for RU B, for an arbitrary input disk u CzV, 
there is an edge between some red dominator of u and some blue dominator of u. 
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Proof. Consider the WVD oi RU B. Without loss of generahty, assume c„ € cell(r) for some r £ R. Now, 
r must be a dominator of u, because r is the closest disk in RU B to Cu- If r does not dominate u, u is not 
dominated by any disk in RU B which contradicts the fact that both R and B arc dominating sets. 

Let b denote a closest blue disk to Cu, that is wvd(c„,6) < wvd(c„,&') for all other blue disks 6'. Note 
that b must dominate u, because if it did not, then no blue disks would dominate u. This would contradict 
the fact that B is a dominating set. Also, note that for any disk d G V such that wvd(c„, d) < wvd(cu, b), 
d must intersect with u. 

If wvd(cti, b) = wvd(c„, r), we are done, since then there is an edge in the dual graph incident on r and 
b. So, let us assume that wvd(c„,&) > wvd(c„,r). 

We will walk from c„ to Cf, along the straight line segment CuCt. The proof strategy is that during this 
walk, we will be crossing red cells and at some point before reaching Cb we will enter a blue cell, in particular, 
cell(6). We must have entered this cell from a red cell cell(r') which shares a boundary with cell(6), and thus 
{r',b} is an edge in our planar graph. Moreover, we will argue that r' necessarily dominates u, completing 
the proof. 

As seen in the proof of Lemma [31 Cf, £ cell(6), and thus we will enter cell(6) at some point in time along 
our walk from c„ to Ch. Let x be the point at which we first enter cell(&). Then x is on the boundary of 
cell(6) and cell(r') for some r' £ RU B. If r' = r, we are done. Otherwise, we have 

wvd(c„, r') < d{cu, x) + wvd(a;, r') ~ d{cu, x) + wvd(a;, b) = wvd(c„, b). 

(Here the strictness of the first inequality comes from our non-degeneracy assumption which implies that 
Cr' cannot lie on the line through Cu and Cf,.) Now, it must be the case that r' G R because wvd(cti,r') < 
wvd(c„, b) and 6 is the closest blue disk to c„. This also implies that r' must dominate u. See Figure[2]for 
an illustration. 




Figure 2: Proof of Lemma [S] The dotted disk is u with center c„ and radius r^. The two red disks r and r' 
are shown as dashed disks with centers Cr and Cr' , respectively. The only blue disk b is shown as a solid disk 
with center Cb- 

Therefore cell(6) and cell(r') share a boundary implying that the edge {&, r'} is in our graph. Moreover, 
b is blue, r' is red, and both dominate w, which completes the proof. □ 

Together, Corollary |4] and Lemma [5] prove Lemma [2] 
2.4 Proof of Theorem [1] 

To show < (1 + e) • we make use of the planar graph separator theorem of Frederickson [11]. This 
argument is similar to the work in [4l[T6] and is only given here for completeness. Given a graph G = (V, i?), 
we denote N{V') for subset of the vertices V to be the set of all vertices in V that share an edge with a 
vertex in V' . 
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Theorem 6 (Frederickson There are constants ci, C2, C3 > 0, such that for any planar graph G = {V, E) 

with n vertices and a parameter r > 1, there is a set X (Z V of size at most cin/^/r, and a partition of 
V \ X into n/r sets Vi, V2, . . . , Vn/r, satisfying: (i) \Vi\ < C2r, (ii) NiVi) C\Vj =0, for i 7^ j , and (Hi) 
\N{Vi)f^X\<c^^. 

We will now show that \B\ < (1 + this is similar to [U [TBI [12]. Let r = b/{c2 + C3) (where b is the 

parameter from the local search algorithm). From Theorem |6|i),(iii), we get \Vi U N{Vi)\ < C2r + c^y/r < b. 
Let Ri = RnVi and Bi = B OVi. Due to the optimality of local search, we must have \Bi\ < \Ri \ + \ N{Vi)\, 
otherwise local search can replace Bi with Ri U N{Vi) to obtain a smaller dominating set, contradicting the 
local optimality of local search. This is why we require that for each d € there is a red dominator of d 
and a blue dominator of d with an edge in the graph. If there were no such edge, then making this swap 
could possibly leave some disks without a dominator. So now we have, 

\B\ < \x\ + J2 m < 1^1 + E i^'i + E ^ i-^i + 



Vb ' 

where c and c' are positive constants. With b a large enough constant times it follows that \B\ < 

il + e)\R\. 



3 The Weighted Dominating Set Case 

In this section, we study a classical generalization of the dominating set problem. Each disk u now has 
an associated rational weight, w„. The goal is to find a dominating set D having the lowest cost, that is, 
wt (D) = J2ueD ""^f small as possible. We will prove the following theorem: 

Theorem 7. Given a disk graph, G = {V, E) of n weighted disks D in the plane, there is a randomized 
algorithm that produces a dominating set V' C V , and wt{V') < 2*^('°s OPT, w.h.p., where OPT denotes 
the cost of an optimal solution. 

The high-level structure of the algorithm is as follows: we first solve a natural linear programming 
relaxation, followed by a randomized rounding step; this step allows us to ignore the weights of the disks 
in the sampling (pruning) stage. In the rounding step, we make several copies of the disks to ensure that 
two properties hold. First, every disk in D is covered by at least n of the copies. Second, the weight of 
the copies is 0{n ■ A*), where A* is the objective function value of an optimal LP solution. Following this 
step, we recursively apply a randomized pruning step where we remove some of the copies according to the 
algorithm given in the proof of Theorem [TU] while ensuring that the remaining copies are a dominating set 
of D. The main goal of the pruning step is to remove some of the copies while approximately preserving 
the ratio of the cost of the remaining copies to the "depth" of the disks in D with respect to the remaining 
copies. We recursively apply the pruning step until the disks in D are covered by only a constant number of 
the remaining copies; the depth of our recursion is 0(log* n). We can then show that the expected weight 
of our final dominating set is at worst 2'^^^°^ - A*. 

First, we define some terms that are used in the remaining part of the section. Given a disk v and a set 
of disks S, we say that v is L-covered by S if there are exactly L disks in S each of which intersects v. In 
other words, neighborhood of v in S has size L. We will make use of the following lemma, which is our main 
contribution to the weighted case: 

Lemma 8. Let S be a set of m disks, and 1 < L < m an integer. Let Q be another (possibly infinite) set of 
disks. There are 0{m- L^) disks of Q that intersect distinct subsets of S each of size at most L. 

Proof. We first define a few concepts that we use in the proof. We focus on subsets S" C S" of size at most 
L and disks of Q whose neighborhood is precisely one of these subsets; let us denote this subset of Q by 



6 



Q' . For a set S" C S* of size at most L, and a pair of disks u,w S Q' , we say that u and w are related if 
they both intersect every disk in S' and no other disk of 5' \ 5", i.e. u n S* = w n 5 = S". So we have an 
equivalence relation on Q' where each equivalence class corresponds to a set S' C S. We wish to bound 
the number of these equivalence classes. Let these subsets of S be {^i, 82, ■ ■ ■ , St}, and correspondingly, t 
equivalence classes {Qi, Q2, ■ ■ ■ , Qt}, where each disk in Qi intersects every disk in Si, and no other disk of 
S \ Si. Consider any set Qi and an arbitrary disk v € Qi- By scaling and/or translating v we can obtain 
a disk v' with the following property: v' has the same neighborhood as all the disks in Qi and is sharing a 
single point with three, two, or one disk in 5*^ and is intersecting all the other disks in more than one point; 
for the cases when v' is touching a single disk in Si, or two disks in Si, we continue to translate and scale v' 
so that it touches two disks outside of Si, or one disk outside of Si, respectively. Without loss of generality, 
we assume that S has four special disks whose borders form the North, South, East, and West boundary, 
respectively, of the region that contains the input disks. We call these special disks TV, S, E, W , respectively. 
Such a transformed disk, v' , that touches exactly three disks is referred to as Vi. We say that a disk d is 
canonical with respect to a set of disks D' if there are three distinct disks in D' such that d intersects the 
three disks at only one point each. Note that each Vi is a canonical disk with respect to the set S. We say 
that a canonical disk v is K-canonical with respect to a set of disks D' if at most k disks from D' intersect 
the interior of v. Therefore, each of the canonical disks Vi that we defined are L-canonical disks. It is easy 
to see that t is within a constant factor of the number of L-canonical disks with respect to S. For each Vi, 
the set of disks that shares exactly one point with it is called the defining set of Vi and every disk of Si that 
shares more than one point with Vi is said to be in the conflict set of Vi. Note that the defining set of Vi has 
at least one disk from Si, but at most two remaining disks can be from outside Si. We will upper bound 
the number of L-canonical disks with respect to S (and hence upper bound t) by choosing a random sample 
S' C S and calculating the expected number of 0-canonical disks with respect to S' . This technique dates 
back to that of Clarkson . 

We pick a random sample S' C S, such that Pr [d £ S'] = j^i while N, S,E,W £ S' with probability 1. 
Let Xd be an indicator variable denoting the event that d is 0-canonical for S' . So, 



Pr[Xd = l]> ( ^ 



1 

1 

L 



1 

1 

L 



> 



1 



where, d is 7-canonical with respect to S and j < L. 

We will show that the maximum number of 0-canonical disks for S' is 0(|5'|). 

Claim 9. For a set S' of disks of size k, the maximum number of 0-canonical disks induced is 0{k). 

Proof. We will bound the number of 0-canonical disks by the number of Voronoi vertices of a weighted 
Voronoi diagram with k sites in which the sites are represented by the k centers of disks in S' , and the 
weight of each site is the radius of the corresponding disk. Every Voronoi vertex is equidistant from the 
disks of the regions sharing that vertex. So each Voronoi vertex in the Voronoi diagram corresponds to the 
center of a disk that touches the boundary of exactly three disks of 5" (disks corresponding to the three 
regions defining that vertex) and does not intersect any other disk of S' . Since the number of Voronoi 
vertices of a Voronoi diagram having k sites is bounded linearly in k, the number of of canonical disks that 
touch three disks of S" are thus bounded linearly in k as well. This leads to the final bound of 0{k) on the 
maximum number of canonical disks that S' admits. □ 

According to the claim, the maximum number of 0-canonical disks for S' is 0(|S"|). So, 



Y,X,<0{\S'\)^^ 



<E[co|5'|] = co-,coGO(1) 



On the other hand, 



E 
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Thus, t < c'mL'^, where c' > cqc G 0(1). 



□ 

We prove the foUowing variant of a theorem of Varadarajan in [19]. 

Theorem 10. Given a disk graph G = {V^E) and set of n weighted disks D C. V in the plane s.t. D 
dominates V , there is a randomized algorithm that produces a subset D' C D, such that for any disk v €z V , 
if V is L-covered in D, then v is at least log L-covered in D' and Pr[d G D'] < £i2Sii. 

Proof. Wc only describe a randomized process that selects a subset, D' of disks such that any disk v E V that 
is covered by D in the range [L, 2L], v is at least log L-covered in D' . Let N,n ~ D, and let Cm denote the 
set of equivalence classes of disks in V such that each class intersects at most 2L disks of D. Note that since 
the disks in one equivalence class of V have the same neighborhood in D, if we obtain a set D' that at least 
logL-covers one disk in that class, then all the disks in that class are also at least log L-covered. Therefore, 
wc can assume we have one representative disk from each class and our goal is to at least log L-cover these 
disks. We use this fact crucially in our analysis. By Lemma [51 \Cm\ < c'-n^L^, — \Nm\- So, there is 
a disk dm that covers at most 2c' classes of Cm- Find such a disk dm S Nm, and recursively compute a 
sequence for Nm-i — Nm \ {dm}, and append the sequence to dm- That is, in the arrangement of Nm-i we 
consider the classes Cm-i whose coverage in Nm-i is at most 2L. The recursion stops when there are fewer 
than L disks remaining, at which point, we compute an arbitrary sequence of the remaining set of disks. 

Let a be the reverse of this sequence, that is, a = [di, ^2, • ■ • , dm)- When considering disk dj, we make 
an instant decision about including it in our cover or not. Call a disk dj £ Nj forced if for some disk v € Cj, 
not including dj will not log L-cover v, whose coverage in Nm is in [L, 2L]. Otherwise, if dj is not forced, we 
add it to D' with probability . We will upper bound the probability of dj being forced - we will show 

that it is at most 0(1/L). 

Observe that if a disk dj is forced because of v, then all the disks dj/ (with j' > j) that cover v are also 
forced, and the number of such disks is at most logL — 1 (otherwise dj won't be forced). So it is sufficient 
to upper bound the probability of a disk di being the first disk forced because of v. Let us denote this event 
by £i{v). Since from among the disks that cover v at most the last logL disks can be forced, the probability 
of one of these log L disks being forced is at most log L times the probability that one of the disks before it 
is the "first" forced disk because of v. We use £i to denote the event that di is the first disk forced because 
of some disk that it covers. 

Consider a disk i; whose coverage is in the range [L, 2L] and di covers it. To bound the probability of 
£i{v) consider the subsequence of cr, called a' , that covers v (so coverage of v is still in [L, 2L]). We will only 
focus on this subsequence since it is the only one which is relevant to £i{v)- Observe that the length of a' 
is in the range [L,2L]. The probability that v forces any disk will then be the sum of the probabilities of 
any of the last logL disks in cr' being a "first" forced disk (by v). Let the prefix of cr' leading up to di be 
called a'i, that is, cr^ = {da^ ,da2, - ■ - , dai=i)- If di is going to be the first forced disk in cr', then i> L — log L. 
Suppose V is the representative disk of one of the 0{L-^) classes of Ci that are covered by di. Since ct-_;^ is 
a subset of disks that cover v, ai-i is also the set of disks that cover v. We will bound the probability that 
an insufficient number of disks were picked from cr-_;^, leading to the forced inclusion of di as the first forced 
disk for v. Let x be the number of disks picked from cr^_j and fix a particular choice of x disks from cr^_x, 
< a; < log L. 



Pr 



X particular i_x — \ 

disks from iy[_-i A disks are 
are picked dropped 



Pr 



Pr 


i—x — 1 


disks arc 




. dropped. 



X particular 
disks from a'^_ 
are picked 

c- logL V / c-logL 
L \ L 
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Since i > L — logL and there are (*^^) ways to choose x disks from 



L-21ogL-l 



< 



Claim 11. 



< (^-^) ( c log£ ] 1 1 cAogL 



e-(L-l) Y ^ ( clog L y , /j^ _ clogL ^ ^ 



: losL \ 2 



_ ( e-e-(L-l)-logL 



,' c-e-(L-l)-logL \ , 4 
max < I } < L 

0<a;<logL [V 1^- X 



Proof. We analyze the continuous function, ^ c e (L^i) iog_L ^ ^ apply the /irst derivative test from ele- 
mentary calculus with respect to a;, x € M. Note that the only value for which the function maximizes 



(L-l}logL'' 



> logL. Restricting the domain, it follows that it maximizes for the boundary point. logL, 



yielding a value of less than L^, when c € 0{1). 
By applying the claim, in the case when x > 0, 

Fr[£,{v)] < L^- 



□ 



< XT, for c> 16 

Note that any disk di/ that occurs before di in a if dii is forced for a disk v' that is not covered by di, 
which forces a disk dk which occurs after di in cr and that dk also covers v, then that event has no bearing 
on the event of di being a first forced disk for v (see Figure |3j). 




Figure 3: di occurs after dii and before d^ in a. di' was forced for w', so d^ is also forced for v' . However, 
dk also covers v. Suppose di becomes forced for v. The fact that dk is already in D' does not increase the 
probability of di's forced inclusion for v. 

So, to upper bound the probability that some dj is a forced disk for a fixed disk v, we sum over all valid 
indices i < j with di being the first forced disk because of v, and obviously there are at most logL of them, 



Pr 



some dj is] — ^ 1 1 



Since there are at most 2c' classes of Cj having coverage in the range [L,2L] that are covered by dj, d 
can be a forced addition for any one of the at most 2c' representative disks. So, 



Pr 



dj is forced 
for some 
disk v£Cj 



< 



2c' 



The probabilistic algorithm finds a dominating set D' C D where the probability of a given disk being in D' 
is at most and each disk v lE V that is covered in the range [L, 2L] by D, is at least logi-covered in 

D' . We repeat the process for points that are between 2L and 4L deep, and so on. Note that the probability 
of a disk being in D' is still the same. □ 
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3.1 Proof of Theorem \7\ 

Let the input instance be a disk graph based on a set of disks D. For any disk d G D, let N[d] denote the 
set of neighbors of d in the graph, inclusive. Consider the following natural LP relaxation for the weighted 
dominating set problem: 

(LP) min ^ WdXd 

deD 

subject to, 

Xd' > 1, Vde D 

d':d'eNld] 

Xd > 0, Vde D 

After solving the LP relaxation, we create a set Dq of disks as follows. For each disk d such that 
Xd > 1/27T., we add Li72;rJ copies of d to Dq. Each copy of d inherits its original cost. For each disk d with 
Xd < l/2n, we don't add any copy to Dq. It is easily verified that wt (-Do) < 2n-X* , where A* is the objective 
function value of the optimal LP solution. Furthermore, we have that each disk d Cz D is n-covered by Dq. 

In the next phase, our algorithm will recursively apply Theorem [TU] to obtain a successively sparse 
dominating set. For the ith application of the theorem, we set Li — logLi_i, for i = 2, 3, . . . , t to obtain a 
set Di C Di-i. For the first application, we set Li = n. 

Fix a disk d and let t be the final application of Theorem [TOl to yield Dt . Let us calculate the expected 
weight of Df. Let Sf denote the event that d & Di. So, 



E[wt(A)] = EPr[^'] 



■ Wd 



(1) 



de-Do 



For a disk d, note that, 



Pr = Pr 



Pr 



t-2 



£ti\/\£t 



Pr 



Pr 



< 



c-log-Lt c-\ogLt-i c-\ogLt-2 c-logL 



t-3 



Lt Lt-i Lt-2 Lt-3 

c- log Lt c- logLt-i c- log Lt-2 e- log-Lf_3 

logLt_i log-Lt_2 log-Lt_3 logit_4 

c*-logLt c*-log(*)Li 



c- log Lj 

U 
c- log -L, 
■ logLi_i 



£i8t\ - Pr [8f\ 
c- logXi 

c-\ogL2 c-log-Li 



logii 



Li Li 

Our recursion stops at a depth t when logLj < 1. Eventually, we get a dominating set -Dj, that is, each disk 
€ y is at least 1-covered in Dt. So, truncating the recursion at a depth of t = log* n, we get that, 

Pr [£t] < f 

Since Li = n and J^deDo "^d — 2't^A*, inequality- ([T]) becomes. 



E [wt [Dt)] < 



Wd < c*2A* 



de-Do 



Using Markov's inequality. 



Pr [wt (Dt) > 2-E[wt (A)] 



1 

< - 
- 2 



Repeating O(logn) independent trials and taking the lowest cost amongst them yields, 

1 



Pr 



> 1 



This proves Theorem [T] 
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Running Time The algorithm requires finding a disk dj in D such that dj intersects 0(m- L^) equivalence 
classes from amongst the subset of disks in V whose neighborhood in D has size at most L. We construct 
these classes, {14, V2, . . . , Vt}, by simply examining a disk v and placing it in Vi such that any disk u £ Vi has 
the exact same set of neighbors in D, and the size of each neighborhood is bounded by L. Since t is bounded 
by Odi?!- L^), a naive approach of comparing the neighborhood of each disk with the neighborhood of each 
of the Vi's, takes 0(|F|- \D\- L'^) time. So naively, the equivalence classes {Vi, V2, . . . , Vt} can be constructed 
in OdVp- \D\- L'^) time. Constructing a set of representative disks U = from the equivalence classes 

takes 0{\D\- L"^) time. Finding a disk dj € D such that dj intersects 0{L'^) disks of U takes 0{\D\'^-L^) 
time. So, constructing the sequence, a, takes 0(|Dp-L'*) time. The quasi-uniform sampling stage takes 
OdD]"^- L^) time for each disk dj because we need to determine if dj is forced for any of the disks in U, any 
one of which may contain some subsequence of disks that occur after dj in a. Since there are a maximum 
of II?! disks that occur after dj, the running time of the sampling stage, given a, is 0{\D\^^■ L^). So, the 
overall running time is 0(|l^p- \D\- + \D\^- L^) for each recursive application of the procedure in the proof 
of Theorem [TUl Since each of is bounded polynomially in n, the number of input disks, and 

given that the depth of the recursion is bounded by 0(log* n) for each independent random experiment, 
and O(logn) bounding the maximum number of random experiments, the over all running time is bounded 
polynomially in n, using a naive approach that does not seek to optimize the running time. 

4 Concluding Remarks and Open Questions 

Given the negative result of Marx |15j which shows that even for the simple case of unweighted unit disk 
graph, an EPTAS for the problem would contradict the exponential time hypothesis [14p|. it is unlikely that 
the dependence of 1/e as an exponent of n on the running time for the PTAS can be improved to, say, 
/(!/£:)• n°(i). However, the running time of the local search PTAS is n*^'^'^ Can this be improved to 
j^O(i/e)? Jn our work, we have made no attempt to improve the running time. 

For the weighted case, we are only able to show a constant integrality gap for the lower bound despite 
numerous attempts. Thus, we believe that the right upper bound for the approximation factor is 0(1). 

Acknowledgments: We thank Saricl Har-Peled and Kasturi Varadarajan for suggesting the use of weighted 
Voronoi diagrams for the unweighted case, and we thank Kasturi Varadarajan for pointing out the connec- 
tion between weighted set cover and weighted dominating set. We also thank Mohammad Salavatipour for 
his support and many valuable discussions. 

References 

[1] C. Ambiihl, T. Erlebach, M. Mihalak, and M. Nunkesser. Constant-factor approximation for minimum- 
weight (connected) dominating sets in unit disk graphs. In APPROX-RANDOM, pages 3-14, 2006. 

[2] Franz Aurenhammer. Voronoi diagrams — a survey of a fundamental geometric data structure. ACM 
Comput. Surv., 23(3):345-405, 1991. 

[3] T. M. Chan. Polynomial-time approximation schemes for packing and piercing fat objects. J. Algorithms, 
46(2):178-189, 2003. 

[4] T. M. Chan and S. Har-Peled. Approximation algorithms for maximum independent set of pseudo-disks. 
In SoCG, pages 333-340, 2009. 

[5] M. Chlebi'k and J. Chlebfkova. Approximation hardness of dominating set problems. In ESA, pages 
192-203, 2004. 

^Marx 1151 actually shows something stronger. 



11 



[6] B. N. Clark, C. J. Colbourn, and D. S. Johnson. Unit disk graphs. Discrete Mathematics, 86(1-3):165- 
177, 1990. 

[7] K. L. Clarkson. Apphcations of random sampling in computational geometry, II. In Symposium on 
Computational Geometry, pages 1-11, 1988. 

[8] T. Erlcbach, K. Janscn, and E. Seidel. Polynomial-time approximation schemes for geometric intersec- 
tion graphs. SIAM J. Comput, 34(6):1302-1323, 2005. 

[9] T. Erlebach and E. J. van Leeuwen. Domination in geometric intersection graphs. In LATIN, pages 
747-758, 2008. 

[10] U. Feige. A threshold of In n for approximating set cover. J. ACM, 45(4):634-652, 1998. 

[11] G. N. Frederickson. Fast algorithms for shortest paths in planar graphs, with applications. SIAM J. 
Comput, 16(6):1004-1022, 1987. 

[12] M. Gibson, G. Kanade, E. Krohn, and K. R. Varadarajan. An approximation scheme for terrain 
guarding. In APPROX-RANDOM, pages 140-148, 2009. 

[13] H. B. Hunt III, M. V. Maratlic, V. Radhakrishnan, S. S. Ravi, D. J. Rosenkrantz, and R. E. Stearns. 
Nc-approximation schemes for np- and pspace-hard problems for geometric graphs. J. Algorithms, 
26(2):238-274, 1998. 

[14] R. Impagliazzo and R. Paturi. On the complexity of k-sat. Comput. Syst. Sci., 62(2):367-375, 2001. 

[15] D. Marx. On the optimality of planar and geometric approximation schemes. In FOCS, pages 338-348, 
2007. 

[16] N. H. Mustafa and S. Ray. PTAS for geometric hitting set problems via local search. In SoCG, pages 
17-22, 2009. 

[17] T. Nieberg, J. Hurink, and W. Kern. Approximation schemes for wireless networks. ACM Transactions 
on Algorithms, 4(4):1~17, 2008. 

[18] S. Pandit, S. Pcmmaraju, and K. Varadarajan. Approximation algorithms for domatic partition. In 
RANDOM-APPROX, pages 312-325, 2009. 

[19] K. Varadarajan. Weighted geometric set cover via quasi-uniform sampling. In STOC, 2010. 

[20] V. V. Vazirani. Approximation Algorithms. Springer- Verlag New York, Inc., New York, NY, USA, 2001. 



12 



